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Abstract. In contrast to recent claims by Enss, Henkel, Picone, and SchoUwock 
[J. Phys. A 37, 10479] it is shown by numerical simulations that the autoresponse 
function of the critical 1+1-dimensional contact process is not in agreement with the 
predictions of local scale invariance. 
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The contact process is a simple lattice model that is often used to describe the spreading 
of an infectious disease [1-3]. It is defined on a hypercubic lattice whose sites can be 
either active or inactive. The model evolves in time by random-sequential updates 
in a way that each active site either spontaneously activates a nearest-neighbor site 
with rate A or it becomes inactive with rate 1. Depending on the control parameter A 
the (infinite) model exhibits a phase transition from a fluctuating active phase into a 
completely inactive absorbing state. This transition is continuous and belongs to the 
universality class of directed percolation (DP). 

Recently Enss, Henkel, Picone, and SchoUwock studied the 1+1-dimensional contact 
process as a generic example of ageing phenomena without detailed balance [4]. In 
particular, they investigated the autoresponse function 



in a critical system starting with a fully occupied lattice at t = 0. The autoresponse 
function quantifies how strongly a local perturbation at position x exerted by an external 
field h{x,s) at a given time s > changes the averaged order parameter (0(5;, s)) at 
some later time t > s. In the contact process the order parameter is the density of active 
sites {ni{t)), where Uiit) = 0, 1 denotes the state of site i at time t, while its conjugated 
field corresponds to spontaneous activation of inactive sites. Hence, the autoresponse 
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function may be determined numerically by simulating a critical contact process starting 
with a fully occupied initial lattice, then activating a random site i at time s (by setting 
ni(s) := 1), and measuring how much this intervention increases the average activity 
{ni{t)) at the same site at a later time t > s. 

It is generally accepted that in scale-invariant systems the response function can 
be expressed in terms of the scaling form 

R{t,s)^s-'-''fR{t/s), (2) 

where is a universal scaling function and the exponent a depends on the scaling 

dimension of the order parameter and the response field. In the case of the 1+1- 
dimensional contact process this exponent is given by a = 25 — 1 ~ —0.681, where 
6 = /3/z/|| is the usual decay exponent of DP (for a definition of DP exponents see e.g. 
Ref. [3]). 

The scaling function /r(^) is usually characterized by certain asymptotic power 
laws. In the case of the contact process these power laws are given by 



MO- , Zl . . (3) 




where \r is the so-called autoresponse exponent and z = denotes the usual 

dynamical critical exponent of DP. The value of the autoresponse exponent was 
estimated numerically by \rI z — 1.76(5). J 

It is important to note that in the critical contact process the two limits of 

(a) first taking s — > oo with ^ = t/s fixed and then sending ^ — > 1, and 

(b) sending s ^ oo while keeping the difference t — s fixed 

give identical results. In the context of ageing phenomena this seems to be surprising 
since in most cases studied so far these limits are known to be different [5]. However, 
if these limits were different in the present case, it would require the existence of a 
crossover from the quasi-stationary to the ageing regime at some typical time scale tc- 
This crossover time tc may be either constant (much larger than the microscopic time 
scale) or it may grow slowly with s such that tc/ s 0. Obviously, in a scale-free system 
such as the contact process at the critical point, where all temporal quantities have to 
scale in the same way, the existence of such a crossover time is impossible, meaning that 
the two limits have to give identical results. In this sense critical systems arc atypical 
examples of ageing phenomena in so far as the crossover from 'equilibrium' to ageing is 
part of the scaling function /ii(^), while in ordinary ageing phenomena it is not. 

I After submission Gambassi [24] proposed that the autoresponse exponent of critical DP may be 
expressed in terms of standard DP exponents as \r/ z = 1 + 5 + d/ z. This sealing relation would give 
the value 1.7921 in one dimension, consistent with the numerical estimate used in the present work. 
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In the case of the critical contact process the equivalence of the two limits implies 
that for ^ ^ 1 the scaling function has to reproduce the behavior of the quasi- 

stationary response function. Since in this limit one measures the response to a local 
activation in the absorbing state, the response function R(t, s) tends to the so-called 
pair-connectedness function c{t, s) (see e.g. [3]) measured at the same positions in space. 
In critical DP this function is known to decay algebraically with the time difference as 
c{t, s) ~ (t — s)"^*^, fixing the exponent in the second line of Eq. (jS)). 

The main point of Ref. [4], on which we will focus in the present work, is to 
apply a generalization of ordinary dynamical scaling to a space-time dependent kind 
of scale invariance, known as the theory of local scale invariance (LSI) [6]. The central 
assumption of this theory is that the linear response function is given by correlation 
functions of so-called quasi-primary fields which transform covariantly under the action 
of a group of local scale transformations. As in the case of conformal invariance, this 
extended symmetry fully determines the functional form of certain universal scaling 
functions. In particular, it is predicted that the autoresponse function is given by the 
explicit expression [6,7] 

(, \ l+a-Xn/z 
-J {t-sr'-'', (4) 

where tq is a normalization constant. This result implies that the scaling function fR{^) 
has to be the form 

fniO oc r'^>(l-l/0-'-% (5) 

which LSI predicts to hold exactly for all ^ > 1. Note that the exponents a and 
\r/ z are not predicted by the theory, they rather depend on the specific model under 
consideration. 

Comparing this prediction with the numerical estimates of /h(0 fo^^ the contact 
process Enns et al. observed a "perfect agreement almost down to t/s = 1" (see 
Fig. 8 of Ref. [4]), concluding that LSI is a symmetry of the critical contact process. 
In a subsequent paper, Ramasco et al. confirmed these findings by Monte-Carlo 
simulations [8]. Moreover, they were able to prove Eq. ^ within a mean field 
approximation of DP, which is expected to be valid in space dimensions d > 4. 

However, there are various reasons to cast doubt on these results. In dimensions 
d < 4 the phase transition of DP is described by an interacting field theory with non- 
trivial loop corrections [9]. To establish LSI as an exact symmetry of DP (analogous 
to conformal invariance in equilibrium models) would be highly non-trivial since the 
field theory would have to be invariant under local scale transformations to all orders of 
perturbation theory. A similar situation is encountered in the two-dimensional kinetic 
Ising model with heat bath dynamics quenched from infinite to critical temperature. 
While Monte-Carlo simulations [7,10] seemed to be in agreement with the predictions of 



Is local scale invariance a generic property of ageing phenomena ? 4 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

1 2 3 4 

t/s 



Figure 1. Data collapse of the autoresponse function R{t,s) in a 1+1- 
dimensional contact process at criticality for s = 1024, 2048, . . . , 16384, using 
the DP exponent 6 ~ 0.15946. The numerical results differ significantly from 
the theoretical prediction of local scale invariance according to Eq. which 
is shown as a dashed line. In order to illustrate how the data complies with 
earlier results, the inset shows the present data on top of a zoomed part of 
Fig. 8 in Ref. [4]. 

LSI, Calabrese and Gambassi [11,12] were able to identify higher-order contributions of 
the corresponding field theory that violate LSI. More recently the existence of corrections 
to LSI in the Ising model could be confirmed numerically [13]. These findings suggest 
to search for a similar failure of this theory in the case of DP. 

Regarding possible violations of local scale invariance, the most sensitive part of 
the predicted scaling function ^ is the knee of the crossover from one power law to 
the other, which takes place at t/s ~ 2. We therefore present a direct measurement 
of R{t, s) in the 1+1-dimensional critical contact process, restricting the simulations to 
the range 1 < t/s < 4. As usual in numerical simulations, scaling can be observed 
only when s, t, and t — s are much larger than the microscopic reference time scale. 
Clearly, such deviations are most pronounced close to t/s = 1 when t — s is small. In 
order to control such errors we generate different data sets, increasing s by powers of 2. 
Collapsing these data sets allows us to get an impression of the range from where on 
scaling sets in. 

These collapsed data sets (generated by a C++ program added as an attachment) 
are shown in Fig. 1. In order to visualize deviations from the prediction of LSI, we 
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plotted R{t,s) s^\l 



1 



{t/ s)^^/^ versus t/s in a double-logarithmic plot. According 



to Eq. (0), LSI predicts this quantity to be constant, as indicated by the horizontal 
dashed line. However, the data collapse clearly demonstrates that there is no such 
coincidence. As can be seen, these deviations are robust and persist as the numerical 
effort is enhanced by increasing s. This means that the predicted scaling form (jH)) 
does not hold in the case of directed percolation, disproving the conjectures of Refs. [4] 
and [8]. 

More recently, Picone, Henkel, and Pleimling proposed a generalized scaling 
form [14, 15] derived from the generators of the Schrodinger algebra. It was argued that 
for applications to ageing, invariance under time translations is actually not needed, 
instead it suffices to require covariance of the response function under the subalgebra 
without time-translations. This assumption then leads to a generalized scaling form 



with an independent exponent a' different from a. This generalized scaling form was 
successfully applied to the one-dimensional Glauber Ising model at zero temperature [14, 
16-18], the OJK approximation [15,19,20], and to the critical Ising spin glass in three 
dimensions [15]. However, in the case of the critical contact process this scaling form is 
inconsistent since for a' 7^ a it does not reproduce the correct asymptotic behavior of 
Eq. © in the limit ^ 1. 

This inconsistency can also be seen numerically. Using the exponents \r/ z and a' 
as fit parameters§ it is possible to find an excellent coincidence with the numerical data 
for a' ~ 0.41 and \r/z ~ 1.76, as demonstrated in the left panel of Fig. 2. However, 
according to the previous argument such a fit has to fail in the limit ^ — 1. As 
demonstrated in the right panel of Fig. 2, this is indeed the case, irrespective of scaling 
corrections due to a microscopic reference time. Similar deviations can be observed in 
LCTMRG calculations [21]. 

The failure of the generalized scaling form Q with a 7^ a' in the contact process 
may be surprising since in various other models it seems to work very well. For 
example, the generalized scaling form was found to hold exactly in the case of the 
one-dimensional Glauber-Ising at zero temperature (and equivalently the annihilation 
process), which is also a scale-free process so that the aforementioned limits are (a) and 
(b) are identical [22,23]. However, in the one- dimensional Ising model at T = the 
magnetic field can only flip spins next to a domain wall whose density is known to decay 
as Therefore, in the limit s 00 with t — s fixed the response function goes to 

zero as -R(t, s) ~ — s)^^/^. Taking the factor into account the generalized 

scaling form can be applied consistently in this case. Contrarily, in the contact process 
the external field can create activity at any inactive site in the bulk so that the response 

§ Note that a cannot be used as a fit parameter because deviations from the DP value 1 + a = 2(5 = 
would immediately destroy collapse of curves for different s. 
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Figure 2. Left: Fit of the generalized scaling function in Eq. on top of 
the collapsed numerical data in the same representation as in Fig. 1. As can 
be seen, there is an excellent coincidence over a wide range of t/s. Right: 
The same data sets plotted against t/s — 1 scaled in such a way that the 
prediction according to Eq. © would lead to a horizontal line. The sharp 
decline for t/s — > 1 illustrates why this scaling form is incompatible with the 
quasi-stationary limit in DP (see text). This decline is stable as s is increased 
so that scaling corrections due to a microscopic reference time as its origin can 
be ruled out. 



function tends to R{t, s) ~ (t — s)"^"^ without an s-dependent prefactor, implying a = a'. 
Similarly, in the 2d or 3d Ising model at criticality, the magnetic field can flip any spin 
in the bulk, hence a numerical flt of the generalized scaling function with a 7^ a' is 
expected to be inconsistent in the limit ^ — 1. 

To summarize, we have shown that the original prediction of LSI in Eq. (0) fails 
in the case of the l+l-dimensional contact process at criticality. In so far the contact 
process is akin to the two-dimensional Ising model at the critical point, where similar 
discrepancies where found explicitly [11-13]. A recently proposed generalized scaling 
form based on LSI without translational invariance produces surprisingly accurate 
results over a wide range of t/s if the exponents are fltted numerically. However, this 
generalized scaling form is inconsistent with the known asymptotic behavior of DP, 
leading to systematic deviations in the limit t/s 1 which cannot be eliminated by 
increasing the numerical effort. These results strongly suggest that LSI in the present 
form (interpreted as a symmetry analogous to conformal invariance in equilibrium which 
determines the scaling function entirely) is not an exact symmetry of l+l-dimensional 
DP at criticality. 

As already conjectured in the context of the 2d Ising model [12, 13], these 
discrepancies indicate that LSI is not a generic property of ageing phenomena, rather this 
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theory seems to be restricted to systems with a Gaussian response, including diffusive 
models [z = 2), off-critical processes, and the mean-field limit (tree level) of certain 
non-equilibrium critical phenomena above the upper critical dimension. 



Note: 

After submission, Enns, Henkel, and Pleimling released a preprint |cond-mat/06052ll^l] 
in which they interpret the deviations shown in the right panel of Fig. 2 as a crossover 
within the scaling function itself that takes place at some value of the quotient C,c = t/s, 
proposing that LSI should hold exactly above C,c- 
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